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Abstract 

The scalar, isovector meson propagator is analyzed in quenched QCD, using 



^ ■ the MQA pole-shifting ansatz to study the chiral limit. In addition to the 

expected short-range exponential falloff characteristic of a heavy scalar meson, 
the propagator also exhibits a longer-range, negative metric contribution which 
becomes pronounced for smaller quark masses. We show that this is a quenched 
chiral loop effect associated with the anomalous structure of the rj propagator 
in quenched QCD. Both the time dependence and the quark mass dependence 
of this effect are well-described by a chiral loop diagram corresponding to an 
rj'-w intermediate state, which is light and effectively of negative norm in the 
quenched approximation. The relevant parameters of the effective Lagrangian 
describing the scalar sector of the quenched theory are determined. 



1 Introduction 



For the foreseeable future, the quenched approximation will continue to play an important 
role in the study of lattice QCD for both practical and theoretical reasons. The practical 
reasons are only too obvious, considering the orders-of-magnitude increase in computer 
resources required for comparable full QCD simulations. But quenched QCD also provides 
a very instructive counterpoint to the full theory which can yield useful insight into the 
effect of quark loops in determining hadron structure. This is particularly true in the light- 
quark limit, where lattice results can be interpreted theoretically with the aid of an effective 
field theory corresponding to quenched chiral perturbation theory. 

The introduction of quenched chiral perturbation theory by Sharpe, Bernard, and 
Goltermann [jl], |2| established the framework for analyzing the chiral behavior of quenched 
QCD. This work focused attention on the particular, anomalous role of the quenched rf 
propagator. The physical rj is approximately a flavor singlet and gets most of its mass 
from the axial anomaly via repeated qq annihilation. In the quenched approximation, a 
single annihilation vertex can appear if both sides of the "hairpin diagram" are attached to 
valence quark lines, but repeated annihilation cannot take place without the closed quark 
loops of the unquenched theory. Instead of cancelling the Goldstone pole in the valence 
quark ("connected") diagram, the single annihilation diagram contributes a term to the rj' 
propagator which has a double Goldstone pole and an overall sign opposite to that of the 
valence diagram. Unlike full QCD where rj loops remain infrared finite in the chiral limit, 
in the quenched theory, the double Goldstone pole of the hairpin term produces additional 
quenched chiral loop (QCL) singularities which alter the chiral behavior of the theory. This 
QCL effect was first observed in lattice results as a deviation from linear behavior of the 
squared pion mass as a function of quark mass H, ||. Recently, the quenched chiral limit 
was extensively investigated in a high-statistics study at j3 = 5.7||. The problem of excep- 
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tional configurations, which had prevented previous studies from probing the smaller quark 
mass region of m n < 350 MeV, was resolved by the pole-shifting ansatz of the modified 
quenched approximation (MQA)H. The QCL effect was observed in both the pseudoscalar 
mass and decay constants at the level expected from a direct study of the rj hairpin mass 
insertion and the topological susceptibility. This established several independent and quan- 
titatively consistent determinations of the quenched chiral log coefficient parameter 5 (or, 
equivalently, of the rf mass insertion 77Jq). 

The spectroscopy of scalar mesons has long been one of the murkier areas of hadron phe- 
nomenology. The well-established scalar mesons are typically at masses of a GeV or more, 
although there have been occasional, fleeting experimental indications of lower mass scalar 
resonances. In the flavor singlet sector, mixing with glueball states further complicates the 
intepretation of experimental data. This is clearly an area where lattice QCD calculations 
can be expected to play a crucial role in the future. High-statistics quenched studies will 
undoubtedly be an important part of this effort. Naively, we might expect that the anoma- 
lous chiral behavior induced by the quenched approximation would be of little importance 
for scalar-meson spectroscopy, since the meson masses involved are all expected to be quite 
heavy and not particularly sensitive to the details of chiral extrapolation. However, as we 
show in this paper, a more careful theoretical analysis combined with a precise numerical 
study of the scalar, isovector meson propagator at very light quark mass reveals perhaps 
the most striking quenched chiral loop effect yet observed. 

The chiral loop effect we discuss here is exhibited very clearly by comparing the scalar, 
isovector propagator for the heaviest and lightest quark masses studied. These are shown 
in Figure |]. For the heaviest quark mass the propagator is a positive and rapidly falling 
exponential, as expected from heavy scalar meson intermediate states. The fit shown is 
a sum of two exponentials with the excited state Oq mass extracted by comparing local- 
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Figure 1: The scalar propagator for the heaviest (« = .1400) and lightest (« = 
.1428) quark masses studied. For heavy quarks (boxes) the propagator is positive 
and exponentially falling. The (upper) solid curve is a two-exponential fit. For light 
quarks (*'s) the propagator exhibits a negative quenched chiral loop effect. 
The lower solid curve is the prediction of one-loop chiral perturbation theory for this term. 
(See Section 2.) 
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local, smeared-local, and smeared-smeared correlators. (See discussion in Section 3.) By 
contrast, the scalar propagator for lighter quarks exhibits a qualitatively different behavior. 
[The lower curve in Figure | is the one-loop quenched xPT prediction for the contribution 
to the scalar propagator of the rf-Tt intermediate state only. See below and Section 2.] This 
unusual behavior has two qualitative features which point clearly to a specific theoretical 
interpretation: (1) The additional component becomes much more prominent for the light- 
est quark masses, and (2) the sign of this component is negative (i.e. opposite to that 
required by positivity of the spectral function). Both of these properties suggest that the 
new component is an effect of the rj'-n intermediate state, which is light in the quenched 
approximation. This contribution arises from the "hairpin + pion" diagram shown in Fig- 
ure |2|(b). This interpretation also explains the fact that this component has a negative 
spectral weight. The hairpin (single mass insertion) diagram appears in the rj' propagator 
with the opposite sign from that of the full propagator. In full QCD the rj'-w intermedi- 
ate state would arise from three different types of quark line diagrams, Figure |2](a) with 
a single vacuum loop going around the entire diagram, Figure |2](b) with a single hairpin 
vertex, and the graphs of type ||(c) with one or more vacuum bubbles inserted into the 
rj' propagator. Of these three, only the hairpin diagram, Figure |2](b), is included in the 
quenched approximation. Since the single hairpin insertion has an overall negative sign 
(corresponding to a positive rj' mass shift), the sign of this contribution is opposite to that 
which would be required by spectral positivity. The prominence of this negative term in our 
measured propagator is a clear example of the unitarity violation induced by the quenched 
approximation. [In the Bernard-Goltermann scalar ghost quark formulation of quenched 
chiral perturbation theory [||, the negative sign of the scalar propagator arises from the 
dominance of a negative metric state consisting of a pair of ghost-quark mesons.] 
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(a) 




(b) 




(c) 



Figure 2: Quark line diagrams which contribute to the rj' — tt intermediate state in 
the scalar isovector propagator. 
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2 Lattice Results for the Scalar Propagator 

In a recent study of quenched chiral logs and the rf propagator || , the MQ A pole-shifting 
ansatz was applied to a set of quark propagators for an ensemble of 300 gauge configurations 
at (3 = 5.7 on a 12 3 x 24 lattice. The fermion action was clover-improved Wilson-Dirac with 
C sw = 1.57. Nine values of quark mass were used, corresponding to a range of hopping 
parameters from n = .1400 to .1428. Valence quark propagators were calculated from both 
local delta-function sources and from exponentially smeared sources in Coulomb gauge. 
(For further details, see Reference ||.) These pole-shifted quark propagators were used to 
calculate the scalar meson propagators considered in this paper. The pion masses obtained 
from this ensemble have values ranging from m^a = 0.245 to 0.603, and are listed in Table 
I. 

It is natural to expect that the scalar iptp correlator in QCD will be dominated by the 
coupling to the lightest scalar meson which is expected to have a mass larger than 1 GeV. 
For the heaviest quark masses studied here, this is in fact the behavior observed for the 
quenched propagator. First, consider the scalar propagator for our heaviest quark mass, 
with m-n-a = .603, corresponding to k = .1400. (From Reference [||, we have k c = 0.14329 for 
clover-improved Wilson-Dirac fermions). The measured values of the local-local, smeared- 
local, and smeared-smeared propagators are shown in the log plot in Figure ||. As expected, 
the plot shows clear evidence of a massive scalar meson. For the smeared-local propagator, 
the time-dependence is reasonably well described for t > 2 by an exponential fit with a 
mass of M s a = 1.25(2). The effective mass plot for the smeared-local propagator is shown 
in Figure || for k = .1400. For comparison, the effective mass for the local- local and 
smeared-smeared propagators are also plotted. Although the smearing function used (an 
exponential with exponent 0.5 in lattice units) was not tuned for this particular problem, it 
does a good job of removing excited state contamination, giving a reasonably flat effective 
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Figure 3: The scalar propagator for the heaviest quark (k = .1400) with local-local 
(*'s), smeared-local (empty boxes), and smeared-smeared (full boxes) sources. 
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Figure 4: Effective mass plot for k = .1400 for local-local (*'s), smeared-local (empty 
boxes), and smeared-smeared (full boxes) correlators. 
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Figure 5: One-loop quenched chiral perturbation theory graph evaluated in Eq. (|2|). 

mass over several time slices from t = 2 to t = 6, and a nearly identical effective mass plot 
for smeared-local and smeared-smeared propagators (indicating an absence of excited states 
in both). 

As pointed out in the Introduction, the behavior of the scalar propagator changes dra- 
matically for the lighter quark mass values. The results for the scalar propagator for our 
lightest quark mass values are shown in Figure ||, with m^a = .245, .267, .307, and .342 cor- 
responding to k values (.1428, .1427, .1425, and .1423) respectively. Here we have plotted 
the local-local ipif) correlators. Instead of exhibiting the behavior expected for scalar meson 
states, these propagators are dominated by a significant negative contribution in the range 
t = 2 to t = 7 which increases for lighter quark masses. We interpret this behavior as a 
clear signal for the negative metric contribution associated with the rj'-ir intermediate state 
discussed in the previous Section. 

The r/-7r loop interpretation of the anomalous component of the scalar propagator is 
made more convincing by noting that the leading chiral behavior of this contribution is 
entirely determined by current algebra in terms of parameters which have already been 
measured with this ensemble in our previous study of the pseudoscalar propagator Q . The 
relevant chiral loop diagram is shown in Figure || with the x denoting the ^'-hairpin mass 
insertion uiq, and the o's representing tptp operator vertices. If we ignore form factor effects 
at these vertices, we can calculate the rj'-ir contribution to the i[)if) propagator in momentum 
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Figure 6: Scalar propagator for k = .1423(«), .1425 (*'s), .1427 (solid squares), and 
.1428 (o). Solid curves are fits to the bubble sum formula, (|l9l). 
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space. The propagator 

A(p) = E^(#W#(0)) (!) 

X 

is given in terms of the one-loop integral B(p), 

A(p) ~ -4rgS(p) (2) 

where 

, . 1 r — > 1 TYIiq , , 

W = yfZ^^ k+ p y + m 2) (jp + m 2)2 ( 3 ) 

Here, 2rg is the t/>V to rf-ir coupling. (In the next section we will discuss including the 
scalar meson intermediate state, along with form- factor /unitarization effects). The pion 
mass values for the hopping parameters used here have already been reported in Reference 
H| and are listed here in Table I. The hairpin insertion mass m§ is also determined in 
Reference ||. For clover improved fermions with C sw = 1.57 the value found in the chiral 
limit (dropping a flavor factor of y/?> included in Table II of that reference) is 

m = .33(2) (4) 

The quark mass dependence of mo was found to be very mild (see Table II of Reference 
H). Here, for simplicity, we will take it to be a mass independent constant given by (|j). 
The only other parameter needed to evaluate the tj'-tv bubble contribution is the coupling 
of the scalar vpip operator to the rj'-n state, 

(0|^|»/7r) (5) 

In the quenched theory in the chiral limit, we can relate the threshold value of this matrix 
element, via a soft-7/ reduction, to the decay constant fp, the vacuum-to-one-pion matrix 
element of the pseudoscalar density. By this argument, we relate the coupling constant in 
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(0) to the m\ vs. m q slope parameter ro determined in Reference ]|], 

r = 1.99(12) (6) 

Thus, in this approximation the effect of the r/-ir bubble, Eq. (|2|), is completely predicted 
by chiral symmetry in terms of the constants (|6|), and the pion masses in Table I. 
The result for the lightest quark mass {k = .1428) is plotted along with the measured 
propagator in Figure [l] (lower curve). Considering that the plotted curve has no adjustable 
parameters, it fits the size and shape of the propagator remarkably well in the time range 
from t = 3 to t = 7. For times t < 3, the measured propagator deviates from the one-loop 
expression in the positive direction, indicating that the spectral function also contains a 
heavy scalar meson. This short-range, positive metric exponential component dominates in 
the heavier quark cases where the chiral loop effect becomes unimportant. The one-loop 
result also explains the pion mass dependence of the anomalous propagator component, as 
shown in Figure |7[ where we have plotted the value of the scalar propagator at various time 
separations as a function of m\. The solid curves are the po Fourier transform at p = of 
the one-bubble integral (|2|), with mo and ro given by (|j) and @. Again the agreement with 
the data is surprisingly good, considering the absence of any adjustable parameters. The 
accurate description of both the time-dependence and the pion mass dependence leaves 
little doubt that, for the lightest quark masses, the scalar propagator over the range of 
times from about t = 3 to 7 is completely dominated by the rj'-ir loop. Furthermore, this 
agreement shows that the soft pion theorem relating the scalar-density matrix element (||) 
to the corresponding pseudoscalar matrix element is well satisfied by the lattice results. 
This is somewhat surprising at (3 = 5.7, where one might have expected substantial chiral 
symmetry violation from finite lattice-spacing effects. 

For short time separations and/or heavier quark masses, an accurate description of the 
scalar propagator requires positive contributions from scalar meson states along with the 
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Figure 7: Pion mass dependence of the scalar propagator for t = 3,4, 5, 6. 
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negative tj'-tt component. In the next Section, we present a derivation of the quenched 
scalar propagator based on a detailed chiral Lagrangian description which will be used to 
fit our data for all quark mass values. The corresponding propagator fits are discussed in 
Section 4. 

3 Quenched Chiral Perturbation Theory 

According to the analysis of Bernard and Golterman, low-energy quenched QCD is described 
by an effective local field theory with degrees of freedom corresponding to meson bound- 
states, including not only ordinary qq mesons, but also q'q, qq', and q'q' mesons, where q' 
is a wrong-statistics ghost quark. Order by order in the rf mass insertion mg, the ghost- 
quark formalism is easily seen to be equivalent to Sharpe's formulation, which begins with a 
U(3) x U(3) chiral Lagrangian combined with rules inferred from the structure of quark-line 
diagrams in the quenched theory. In this paper, as in our previous analysis Q we follow the 
latter approach, invoking a chiral Lagrangian having valence pion and rj' degrees of freedom. 
To describe quenched QCD, the chiral Lagrangian is supplemented with rules which reflect 
the suppression of internal quark loops. As a result, all normal pion loops are suppressed 
except for cactus diagrams involving the r/ meson connected by single insertions of the 
hairpin mass term. Using this procedure we obtained a consistent fit to our lattice data 
for the pseudoscalar valence and r( hairpin propagators and determined the relevant chiral 
Lagrangian parameters. We now extend this analysis to the isovector, scalar propagator. 

The appropriate effective field theory must now include a multiplet of scalar mesons in 
addition to the pion and if degrees of freedom. Using standard chiral Lagrangian methods, 
the heavy scalar mesons are described by scalar fields transforming nonlinearly under chiral 
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symmetry rotations. The resulting chiral Lagrangian is 

f 2 f 2 
C = J —tr{dUdU ] }+ J —tr{ X ] U + U^x} 

+ jtr{DaDa} - ^m 2 s tr{aa} + g s tr{a^UdU ] dU \/U ] } 

+ fsHxWUvVU + xVU^VU^ + C^irpin (7) 

where I? is a chiral covariant derivative, m s the sigma mass, g s the strength of the chiral 
invariant coupling of sigma mesons to pions and f s the strength of the corresponding chiral 
symmetry breaking interactions. Here we denote the scalar field by a. (Elsewhere, when 
referring to particle states, we use the conventional ao and ag to denote ground-state and 
excited-state mesons.) The chiral field matrix U contains the pion and rf degrees of freedom. 
The rj mass term is given by the hairpin Lagrangian, 

£hair P in = - ^(/ 2 /8) [itr ln(C^t) _ itrln(U)} 2 . (8) 

We have not included interactions involving the couplings L5 and Lg used in our previous 
analysis of higher order terms in the sigma potential as they will have only small effects on 
our analysis of the isovector scalar propagator. 

Scalar and pseudoscalar quark densities are represented by meson operators in the ef- 
fective field theory and can be determined from the dependence of the chiral Lagrangian on 
the spurion field x- Hence the isovector scalar density operator is given by 

^2*1 = -\r f 2 (U + U^ 12 

-2r f a (y/U<Ty/U + v^W^ t )i2 (9) 

with 

x = x f = ml = 2r m quark . (10) 
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Our previous fits to the pseudoscalar propagator determined the following values for the 
chiral Lagrangian parameters, 

/ = 0.1066(24) 
m Q = 0.33(2) 

r = 1.99(12). (11) 

To lowest order, the local isovector scalar density may be expanded in terms of the pion, 
rf and sigma fields, 

¥i = 2r (7T+7/) " 4v / 2r / s (a+). (12) 
At tree level, only the sigma propagator gives a positive contribution to the scalar propagator 

<^*! W& 2 >= 32r% ffPr. (13) 

where P a is the ao meson propagator. The rj'-n term can only contribute via a meson loop. 
However, the rf propagator in the loop is the rf hairpin propagator with a single insertion 
of the hairpin mass, m§ (see Figure §). Including only these terms, the scalar propagator is 

< W 2 y 1 W[V 2 >= 32^/ s 2 P CT + 4rgS hp (14) 

Here Bhp is the hairpin bubble 

B hp = FT {yfY< ((A . + p} 2 + m 2) ( fc 2 + ^2 )2 } ( 15 ) 

and FT means the Fourier transform on p. As discussed in the Introduction, the negative 
sign of the r( propagator in loops is dictated by the quenching and implies that this meson 
loop makes a negative contribution to the scalar propagator. It should also be noted that 
the rf propagator is more infrared singular than the normal valence pion propagator. This 
has a dramatic effect on the scalar propagator. In the light quark limit, the quenched theory 
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predicts that large negative chiral loop effects should dominate the scalar propagator. Our 
lattice data give convincing evidence for these effects. 

A complete analysis of the scalar propagator must go beyond the lowest order terms 
given above. The chiral Lagrangian predicts couplings between the scalar meson and rf-ir 
states as well as rf-n rescattering interactions generated by the symmetry breaking terms 
that give mass to the pion. These higher order terms involve rf-ir bubble diagrams and 
can be resummed into a closed form. Since the chiral invariant interactions between the a 
and rf-ir states include derivative couplings, bubbles with derivative couplings at the arfir 
vertex must, in principle, also be included. These terms are less infrared singular than 
the nonderivative terms and are not expected to play a major role in explaining the large 
quenched effects seen in the data. In the following analysis, we will neglect such derivative 
coupling terms. The formulas and results given below are for the equal quark mass case; 
generalizations to the unequal mass case are straightforward. 

As is always the case in propagator studies, the presence of excited states in the spectral 
function complicates the analysis. As we discussed in Section II, in addition to the propaga- 
tor for local ipip operators, we have calculated propagators with a smeared source at one or 
both ends. In our previous analysis of the pseudoscalar channel |5| , a full multi-state fitting 
procedure was employed. Here we take a somewhat simpler approach which is suggested 
by features of the data. We model the spectral function in terms of three components: (1) 
the rf-ir state, (2) a ground-state scalar meson ao, and (3) an excited-state scalar meson 
<2q. We perform an analysis of the local-local propagator by first using information from 
the smeared-local and smeared-smeared propagators to remove the excited-state Oq scalar 
meson component and then fitting to a formula obtained by resumming all repeated bubble 
graphs involving the rf-ir state and the ground-state ao scalar meson, as depicted in Figure 
|8[ To estimate and remove the contribution of the excited <Zq from the local propagator, 



18 



we consider first the case with the heaviest quark, k = .1400. Within statistics, all of the 
scalar propagators for smeared and local sources exhibit the same time-dependence beyond 
t = 3. In fact, the smeared-local and smeared-smeared propagators are equal, up to an 
overall factor, at all time slices, as evidenced by the effective mass plot in Figure ||. For the 
heavier quark mass values, the tj'-tt loop is relatively unimportant, and the smeared and 
local propagators can be analyzed in terms of a ground-state and excited-state scalar meson. 
The fact that the smeared-local and smeared-smeared propagators, Asx(t) and Ags(t) are 
proportional and that both give flat effective mass plots indicates that neither propagator 
has a significant excited-state component. On the other hand, the local-local propagator 
clearly has a significant Oq component, as seen in the effective mass plot, Figure |j. 

It is straightforward to estimate the Oq contribution by using the (essentially time- 
independent) ratio Asl/A-ss to determine the excited state component of All,. 

A L L{t) - ^-^Asx(t) ~ C ao * exp(-m ao *f) (16) 
Ass 

Numerically, the above procedure applied to the k = .1400 propagators gives excited state 
parameters 

C a * = 0.275(30) , m al = 1.86(10) (17) 

The large mass obtained emphasizes that the excited state contribution is only important 
at very short time separations. Similar results are obtained for k = .1405 and .1410, with 
no significant observed dependence on quark mass. With this information, we can separate 
off the excited state component of the local-local propagator: 

A LL (t) = C a * exp(-rn a *t) + A LL (t) (18) 

where All{P) includes contributions from the rj'-n state and from the ground state scalar 
meson, and will be modelled by the bubble sum depicted in Figure For the lighter quark 
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masses, the tj'-tt state becomes important, and it is more difficult to disentangle the ao and 
ag contributions. In the fits to the local-local propagator described in the next section, 
we have neglected the quark mass dependence of the excited state term and used the 
estimate fllTp for all quark masses. 

By resumming the multiple bubble and scalar propagator graphs shown in Figure [8|, the 
full propagator may be written in terms of the one- loop bubble function B(p) defined in 
equation (g), 

A(p) = -4r 2 B ren + 32r 2 (/ s - f s (2m 2 Jf 2 )B ren ) (/. - f s (2m 2 Jf 2 )B r£n ) P a0 (19) 

where f s is the coupling of the scalar meson to the ipip density, B ren (p) is the resummed 
multiple-bubble function, 

Tjren B(p) (nri\ 

~ 1 + {2ml/P)B {2U) 

and P a Q is the resummed ao propagator, 

Pao =P 2 + m 2 s + 32(f s m 2 Jf 2 ) 2 B ren (21) 

Note that the negative sign associated with the hairpin insertion has been taken outside of 
the bubble function, i.e. the functions B(p) and B ren {p) are both positive. 

We use the bubble function B(p) defined for the lattice version of the full chiral La- 
grangian with the same lattice spacing as the quark lattice. In (20,21), the lattice four- 
momentum is defined by p = 2sin(p/2) as appropriate for bosonic propagators. If we 
are to use the pion masses in Table I, which were determined from fitting the t-space 
pion propagators, then the pion masses used in lattice propagators above are given by 
rh n = 2 sinh(m vr /2). 

This analysis includes all contributions to the scalar propagator which can be computed 
by direct rj'-ir bubble summation and mixing with the ao scalar meson state. We use these 
expressions in our analysis of the lattice data on the isovector, scalar propagator. 
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The quenched theory does introduce additional contributions which could modify the 
behavior of the bubble terms. For example, pions can interact through double hairpin 
exchange diagrams whose infrared behavior is doubly enhanced. Therefore, virtual rf pro- 
cesses do not decouple in the infrared limit as is normally expected for soft pion interactions 
but can generate long range forces between mesons in the quenched theory. Preliminary 
estimates of these effects indicate they are negligible for the range of pion masses and lattice 
volume considered in this paper and they have not been included in our analysis. 



4 Global Fits and Determination of Chiral La- 
grangian Parameters 

In Section 2 we showed that the data for the scalar propagator for heavier quark masses 
was dominated by the scalar meson ao intermediate state (after removing the excited ag 
contribution, as discussed in that Section), while for light masses, it is dominated by the rj'-ir 
intermediate state. To get a consistent fit for all masses, we have found that the resummed 
formula fll9D (Figure ||) produces the most stable results. (The simpler alternative of using 
the sum of a scalar meson pole plus the single bubble function gives fits with only slightly 
worse x 2 ' s but the resulting fit parameters f s and m s are less stable as a function of 
quark mass.) All of the fits discussed in this section are based on the formula (19) with the 



parameters ro and / fixed at their previously determined values (|il|). We present the results 
of two different fitting procedures. In the first set of fits (Table I), we take the hairpin mass 
insertion parameter rriQ to be fixed at the value ([ll]) and use the scalar propagator fits to 
extract the scalar meson Lagrangian parameters f s and m s . In the second set of fits (Table 
II), we let mo be a fit paramter and investigate how well the scalar propagator data by itself 
determines the value of the hairpin insertion. 
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Table 1: Fit parameters for the scalar propagator, with ttlq fixed. The fit range is 
f=[l-6]. 



K 


trip 


/. 


m s 


X 2 /dof 


.1400 


.603(2) 


.0544(10) 


1.20(2) 


1.7/4 


.1405 


.556(2) 


.0548(10) 


1.19(2) 


2.0/4 


.1410 


.505(2) 


.0550(11) 


1.17(2) 


1.9/4 


.1415 


.450(3) 


.0554(12) 


1.16(3) 


1.7/4 


.1420 


.386(3) 


.0559(14) 


1.16(4) 


1.7/4 


.1423 


.342(4) 


.0565(19) 


1.18(5) 


2.3/4 


.1425 


.307(4) 


.0569(25) 


1.19(7) 


3.7/4 


.1427 


.267(5) 


.0568(38) 


1.20(12) 


5.6/4 


.1428 


.245(6) 


.0571(54) 


1.29(18) 


6.2/4 



The fits which extract the scalar parameters, with uiq held fixed, are given in Table |l|. 
We see that the value of the ground state scalar mass parameter m s is rather well determined 
for the heaviest quark masses, but becomes less accurate for lighter masses, where the rf-n 
state dominates. A linear fit to the scalar mass values in Table |l| gives a value of 

m s = 1.14(7) (22) 

in the chiral limit. Using the charmonium scale a -1 = 1.18 GeV, this gives a scalar meson 
mass of 1.34(9) GeV. The scalar decay constant f s in the chiral limit is 

fs = -057(3) (23) 

In addition to extracting scalar parameters, it is interesting to carry out fits to the scalar 
propagator with rriQ as a fit parameter. This procedure emphasizes the fact that the scalar 
propagator data alone provides a fairly accurate estimate of mo which is independent of, 
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Table 2: Fit parameters for the scalar propagator, with m s fixed to values shown. 
The fit range is t — [1 — 6]. 
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m s 


/, 
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X 2 /dof 


.1400 


1.200 


.0568(12) 


.336(22) 


2.1/4 


.1405 


1.185 


.0570(12) 


.335(25) 


2.3/4 


.1410 


1.170 


.0571(13) 


.335(27) 


2.3/4 


.1415 


1.156 


.0574(13) 


.338(30) 


2.1/4 


.1420 


1.142 


.0575(14) 


.346(33) 


2.3/4 


.1423 


1.133 


.0578(15) 


.358(36) 


3.0/4 


.1425 


1.128 


.0576(16) 


.360(43) 


4.5/4 


.1427 


1.122 


.0564(18) 


.332(51) 


6.4/4 


.1428 


1.119 


.0555(20) 


.298(64) 


6.7/4 



and consistent with, the previous estimates in Reference ||. As we saw in the previous fits, 
the scalar mass m s is well-determined for the heaviest quarks, but is poorly determined at 
the light-quark end. To obtain stable fits with uiq as a fit parameter, we fix m s to be given 
by a linear fit obtained from the three heaviest quark values, m s = 1.106 + 1.145m g . where 
m q is the bare quark mass (k" 1 — k~ 1 )/2 in lattice units. In this set of fits, f s is found to 



be approximately independent of quark mass and is well-approximated by (23) for all quark 
masses. The values of the rf mass parameter mo obtained from these fits are given in Table 
|2] and plotted in Figure |9[ A linear fit to this data gives 



mo = 0.34(4) 



(24) 



in the chiral limit. This result, obtained from the scalar propagator data alone, is quite 
consistent with our previous estimates which included direct r( hairpin measurements, a 
topological susceptibility calculation, and QCL effects in the pseudoscalar channel 0. 
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Figure 9: Results for the rf hairpin mass insertion extracted from fits to the scalar 
propagator. The solid line is the value obtained from direct analysis of the hairpin 
diagram in Reference || (after removing a \[Z flavor factor included in the definition 
of mo in that reference). 
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One feature of the scalar propagator which appears to be somewhat inconsistent with 
our theoretical model is the behavior at larger time separations, t > 7. For the lightest 
quarks, we saw that the propagator was well-described in the region t = 3 to 6 by the one- 
loop rj'-TT bubble. However, beyond t = 6, the propagator appears to go to zero (from below) 
more rapidly than the theoretical curve predicted by either the one-loop calculation or the 



bubble-sum formula, as shown for k = .1427 in Figure 10. If it is a significant departure 
from the chiral symmetry prediction, it could be an indication of strong rj'-ir interactions 
near threshold. Such interactions would normally be ruled out by soft-pion arguments, but 
may arise due either to explicit chiral symmetry breaking of the Wilson-Dirac action or 
because of higher order quenched effects. For the heaviest quarks, we also find that the 
propagator appears to have an additional positive contribution for t > 7. For example, we 
can obtain a good fit (to both the heavy and light mass data at large time) by adding a 
small 7T propagator term to the fit model. However, the origin of such a contribution is not 
presently understood within our model. 



5 Discussion 

We have presented convincing evidence that the quenched scalar, isovector meson propaga- 
tor exhibits a prominent "quenched chiral loop" (QCL) artifact at light quark mass which 
is properly interpreted as the effect of an tj'-tt intermediate state. QCL effects that have 
been observed in other quantities ]5|, |[] have been rather subtle, involving relatively small 
deviations from full QCD chiral behavior. By contrast, the QCL effect discussed here is an 
unmistakable chiral power singularity, and in fact dominates the scalar correlator for the 
lightest quark masses we have studied. The negative spectral weight of this contribution is 
a particularly vivid illustration of the oft-stated truism that the quenched approximation 
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Figure 10: Comparison of scalar propagator for k = .1427 with the bubble sum 
formula (|i~9"D fitted to the interval t=l-6. 
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violates unitarity. 

The results presented here also provide further evidence of the effectiveness of the pole- 
shifting ansatz of the modified quenched approximation. It is only after pole-shifting that 
we are able to study the lightest quark masses for which the rj'-ir state becomes obvious. 
It should be remarked, however, that even for the heavier quark masses, the QCL effect 
is not negligible, and must be included in the fitting function in order to draw correct 
conclusions about the scalar meson mass. For example, fitting the scalar propagator at the 
heavier quark masses with a simple exponential and ignoring the rj'-ir contribution would 
lead to the erroneous conclusion that the scalar meson becomes heavier as the quark mass 
gets lighter. The resulting effect on the chiral extrapolation is at least as severe as typical 
QCL effects in the quenched light hadron spectrum. For pion masses below about 400 
MeV, the rj'-TT contribution becomes dominant, making the scalar meson mass difficult to 
determine accurately. Nevertheless, by combining lattice calculations with quenched chiral 
perturbation theory, more extensive studies of scalar meson spectroscopy in the quenched 
approximation should be feasible. Finally, the prominence of the QCL effect in the scalar 
valence propagator suggests that this propagator be added to the list of standard bellwether 
quantities (such as topological susceptibility, r/ mass, and long-range static potential) which 
are expected to be particularly sensitive to the difference between quenched and full QCD 
gauge ensembles. 
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